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ABSTRACT 


This document discusses some of the implications of the Kirchhoff- 
Kottler formulation for discontinuous surface distributions. The 
intrinsic impedance of free space is recovered from the equations for 
the electric and magnetic fields, and the gain dependence of the antenna 
on frequency and geometry is identified. The frequency dependence of 
each integral is tabulated, mid each is associated with a physical quantity 
on the reflecting surface. It is shown that the Kottler contour integral 
vanishes for closed surfaces, independent of field values on the boundary. 
The behavior of Vi//, the gradient of the solution to the wave equation, is 
examined in a local context on the antenna surface to illustrate that 
charge distribution can give rise to transverse field components in the 
intermediate near-field region. Finally, derivations originally due to 
(1) J. A. Stratton and (2) M. I. Sancer are reviewed to bring out the fact 
that the Kirchhoff-Kottler formulation inherently satisfies the radiation 
conditions and vanishes the radial fields of scattering surfaces in the 
limit as the observer recedes to infinity. This important conclusion is 
established by two unique approaches. The development of Sancer»also 
demonstrates that the Kottler contour integral can be evolved without 
any appeal to physical intuition, for the case of discontinuous illumination 
distributions, via the dyadic Green’s function. 
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THE CALCULATION OF ELECTROMAGNETIC FIELDS 
BY THE KIRCHHOFF-KOTT LER METHOD 


INTRODUCTION 

The starting-point of this document is the Kirchhoff-Kottler formulation for 
computing backscattered electric and mrgnetic fields from surfaces carrying 
charge and current distributions. Although standard textbooks 1 have carried 
this formulation for many years, the analytical difficulties associated with the 
evaluation of the integrals have delayed application to many practical antenna 
configurations. The advent of high-s|>eed digital electronic computers has made 
possible the numerical evaluation of the otherwise intractable integrals, and 
has brought about renewed interest and a re-examination of the theory. 

A few' dates are interesting and significant. Maxwell's equations appeared around 
1865, and Kirchhoff's theory was presented about 1882. By 1923 Kottler 2 had 
annexed his contour integral which, when taken in conjunction with Kirchhoff's 
results formed a theory which satisfied Maxwell's equations. More up-to-date 
discussion of the Kirchhoff-Kottler formulation appears in Physical Review by 
Stratton and Chu' in 1939, the textbook by Stratton 1 in 1941, and the well-known 
microw'ave antenna text bv Silver 1 in 1949. Jones 4 mentions the Kottler in- 
tegral again in the representation of the electromagnetic field in his 19G4 text. 

Up to this time, it appears that the Kirchhoff integrals over current and charge 
distributions are regarded as basic, and the Kottler contour integral on boundary- 
line charge is considered a "correction" term or approximation. A 1964 textbook 
by Van Bladel/ which makes extensive use of the dyadic Green's function in the 
calculation of fields, comments on the linear charge density, and provides back- 
ground for a unified approach to the diffraction formulation. In Radio Science, 
1968, Sancer 6 presents a derivation of the general vector Kirchhoff equations, 
via the free-space dyadic Green's function, which intrinsically includes Kottler's 
boundary-line charge for open surfaces. The latter derivation is in a sense com- 
plete, and does not appeal to separate physical argument or intuition in achieving 
its objectives. 

The Kirchhoff-Kottler formulation is usually presented in the notation of vector 
analysis, and constitutes a rather formidable array of symbols. It is abstract in 


^Ref. 1, p. 
2 Ref. 3 
*Ref. 4 
'Ref. 5, p. 
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6 Ref. 7 


469, Ref. 2, p. 160 
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that the general formulation simply postulates existing fields (Ej, H,) which are 
to be integrated over some unspecified surface. The computer programmer 
usually notes that comple x vector s are used throughout so that each of the eight 
original integrals decomposes into six integrals in the real domain (appropriately 
sorted according to their Cartesian vector components and the real-imaginary 
quadrature). Sometimes the program language and computer system allow direct 
vector operations and complex algebra. Nowhere, in the original formulation, 

Is there an explicit connective between the sources and the illuminated surface. 1 
The boundary conditions at the surface are not specified in the general formula- 
tion. In short, a theory descriptive of a physical process lias been made avail- 
able, not a flexible empirical formula which is adapted to a large class of prob- 
lems on the basis of existing knowledge of antenna systems as is sometimes 
suggested. The theory is used exactly as it stands, without modification of any 
kind, and the test of its validity lies in comparison with other analyses and care- 
ful measurements on actual physical systems. 

One way to remove the abstract character of the general formulation is to begin 
to specialize its application and examine the simplified result. I’robably the most 
natural and practical assumption is that of a perfect conductor, leading to the 
concept of a sheet current and a charge distribution on the surface of the idealized 
conductor. 2 The assumption of an infinitely remote observer affords another 
opportunity to gain familiarity with the notation and the Kirchhoff-Kottler formula- 
tion. When these assumptions are made the gain dependence on frequency and 
geometry can be seen clearly for both E( x ' , y ' , z ' ) and H( x ' . y ' , z * ) , An 
identification of the radial and transverse spherical components of the composite 
field is easily made for the case of the distant observer. 

Returning to the general formulation, and assuming a perfectly conducting sur 
fa ?e, the intermediate near-field can be studied. The function of the term V0, 
where 4 > ~ c~* kr /r is taken as a solution to the wave equation, in a local context 
on the antenna surface can be related to an intuitive concept of fields due to 
point-charge distributions. Finally the question of whether or not the complete 
Kirchhoff-Kottler formulation satisfies the radiation condition 3 can be answered 
without resorting to numerical evaluation or reference to specific illuminations 
or reflector geometries. The vanishing of the Kottler integral for closed surfaces, 
under arbitrary illumination conditions, is obtained simultaneously while resolving 
the preceding question. It has been found that these exercises are invaluable in 
purging errors from a freshly written computer program and contribute to the 
overall appreciation of the Kirchhoff-Kottler diffraction theory. Most striking, 
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perhaps, is the intern:) 1 self-consistency and close a<lhorence to well-established 
physical principles that emerge through application of the theory. 


GENERAL FORMULATION 

The backscattered electric and magnetic fields of the general Ki rchhof i -Kottler 
formulation are: 


E(x\ y\ z') = - V0Hj 


i\Z 


© 

(Kott lei ) 


“ 4 ^ J H,) 4 > * (n xEj) xV*//+ (ii • E,) Vi//] ds 


® ® 

( Ki rchhof f ) 




H( x* , y\ z') = 


j ix 4n | 
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(Kottlcr) 


+ l J [ j* e (n x E,) 0 - (n * H, ) x\ty - ( n • Hj ds 


® ® 

(Ki rchhof f ) 
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whore 


0 




In the above, E, and II, are electric and magnetic distributions which exist in 
space and may or may not be associated with an actual physica) urface. When 
a material surface is chosen, appropriate boundary conditions are applied to 
evaluate E,, H, in terms of the incident fields E |t ll ( . Vectors such as n, E,, H, , 
t\l t and V0 are ordinarily s|>ecified by their Cartesian components and associated 
with a definite point of application in the pro un. They are bound, 2 not to the 
origin, but to the surface and are taken in a local context. These vectors are 
part of an integration on the four-dimensional manifold (l H , l y , l t , j) where the 
Cartesian basis vectors are fixed in space. The resulting vector fields E(x' , y', /' ), 
H(x' , y ' , z' ) are usually rewritten in terms of projections on the moving spherical 
triad (l r , 1^. 1,) and a phase angle. 

SPECIALIZATION OF THE GENERAL FORMULATION 

When the general formulation is applied to the class of problems for which the 
illumination distributions reside on perfectly conducting surfaces, the Kirchhoff- 
Kottler theory reduces to: 3 


E(x\ 
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^Silver’s conventions are used for \pin the above* 

"Ref. 8, pp. 33-34, Ref. 9, pp. V6. These references distinguish several kinds of vectors. 

^The component fields are radial and transverse, as indicated, for an infinitely remote observer. 
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as before. 

The simplification is due to the boundary conditions' 
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Physically, the antenna surface is considered locally plane, and the normal com- 
ponents of H, and the tangential components of E, vanish at the interface. In 
effect, the normal E, and tangential H. fields are converted into electric charge 
distributions, and electric sheet currents flowing along the local tangent to the 
surface. The contour integral can be regarded as a discontinuity of the surface 
sheet current resulting in a surfeit or deficit of electric charge. 


The vector mechanics by which the Kirchhoff-Kottler formulation achieves 
resu’ts satisfying physical intuition about the problem may be worth reviewing. 
Vector dot- and cross-products can serve as mnemonic devices for identifying 
sheet currents lying in the reflecting surface, surface charges, and boundary 
charges. It is fundamental that a dot-product amounts to a projection between 
two vectors and a cross-product always results in a vector w'hich is orthogonal 
to the given pair of vectors. From this it follows that (n • E,) Vi// of integral @ 
is due to that part of E, w'hich is orthogonal to the surface and induces charge. 
Since ( n • Ej ) V<// is radial, this component of the total field acts radially on an 


*Ref. I, pp. 34-3*' 
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observer, as do all fields from charges. 1 Also the term (iixH,) must be a 
vector which is orthogonal to both n and H,. If iixH, is orthogonal to n it 
must be in the tangent plane of the surface, and is suited for describing a sheet 
current due to tangential H, as in integral ( 2 ). Depending on the orientation of n 
over the surface, ( n x H, )can obviously exhibit either radial ji transverse fields 
at the point of the observer. In integral (X the sheet current is crossed, vcctor- 
ially, into V</>, which is purely radial. Them the integrand (fix ii 1 )xV<// is, 
necessarily purely transverse, as is intuiti vely required of the scattered magnetic 
field. Finally, the contribution of integral CO * s radial since H, • d~ is a scalar 
and Vi p is purely radial. Since d t lies in the surface along the contour and is a 
tangent to it in a local context, the magnitude of H, • dT depends on the relative 
aspect between H, and d i according to the projection rule. Obviously the tangential 
part of H, affects the magnitude of the charge on the contour as stated previously. 

The preceding conclusions about the radial and transverse nature of the field 
components were obtained on the assumption that the observer was infinitely 
remote and the surface was perfectly conducting. It could be seen that H(x\ y'. z') 
was transverse since any vector crossed into a radial Vi// is transverse, and 
only a single integral provided the solution for the scattered magnetic field. The 
scattered electric field E(x' , y' . z' ) contained one integral with both transverse 
and radial components, and two other integrals were shown to contribute only 
radial fields. It is not a priori evident that the three radial components combine 
cooperatively, in general, in the Kirchhoff-Kottler theory and annihilate to satisfy 
the radiation condition. 2 On the assumption that the radial fields might vanish, 
the remaining transverse fields should yield a w-ell-knowm ratio, the intrinsic 
impedance oi free space. 


iEi'x', y', »‘)| 
!R(x\ y\ i')l 




(n x H) x V0 ds 



(nx Hj) 'p ds 
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,s i 


-iu f 

(nx Hj) x 1 r 0ds 
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( f ree space ) 


Ref. 1, Appendix I, p. 602 gives the dimensions of the physical quantities that also assist in 
identifying the behavior of the integrands. The physical constants go, /jl, and € are not included 
in the discussion above, and only the vector character of the integrands is examined. 

“The discussion of the proof is deferred to the latter sections of this document. See Appendix A 
and Appendix B. 
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The cross-product involving l r in the denominator causes no difficulty here. 
Since it has been assumed, tentatively, that the fields are purely transverse for 
the infinitely remote observer, it follows that 



From this preliminary analysis it appears that the Kirchhoff-Kottler theory can 
provide a diverging spherical backscattered wave, subject to additional analysis 
concerning the vanishing of the radial fields. 

Additional confidence in the description of the backscattered fields can be ob- 
tained by postulating a pair of incident fields E ( , H i which give a Poynting vector 
P t = E t x ii i tha* is directed to ward some scatterer. A locally plane segment of 
the scatterer can be taken, and the directions of the scattered fields are assigned 
by the integrals © and © of the theory, with due regard for signs. It is easy to 
verify that the scattered fields give a Poynting vector P s = E(x', y', z') x H(x', y\ z') 
away from the reflector. This result, although somewhat trivial in appearance, 
is important since subsequent applications utilize the backscattered fields to 
illuminate other surfaces (i.e. dual reflector systems such as those of Cassegrain 
and Gregory). 1 The E(x' , y' , z'),H(x', y', z') comprise a ne\v illumination 
function, effectively, and determine a new E t , H, , but the incident fields are not, 
in general, related simply to one another as in the case of a spherical wave. 


•Ref. 1, pp. 131-137 
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A subject of considerable interest in the antenna field is directive gain, which is 
defined as 1 


G(<9, 0) 


if r 

J 0"O J <t> m 0 


P(0. 0) 


P(&, 0 ) sin & d 0 J 0 


P, («,«) + P 4 ,( 6 ,*) 


iff 

•'o J o 


P* ( 0 , <p) ds + 


±JT 

•'o •'o 


( 0 , 0 ) ds 


This topic will not be discussed at length but it is interesting to discover how 
gain variations are introduced in the Kirchhoff-Kottler formulation. 

A special class of scatterers is considered, initially, for which all the differential 
contributions of integrals (5) or (?) arrive in-phase at the point of the observer. 
Focal-point fed paraboloids, observed on-axis, are familiar examples. The 
directive gain of such reflectors is usually given by 


G (6, 0) 



where c is an efficiency factor related to illumination edge taper. For a fixed 
geometry diameter d is fixed and 


G(0, 0) * U? 


In the Kirchhoff-Kottler formulation, considering both the electric and magnetic 
fields, the same frequency dependence is required by intuition. 


] Ref. 9, pp. 179-180, Ref. 2, p. 90 
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Since a cophased aperture exists by assumption, 


G(<9, <t>) ? | E( x' , y ' , z ' ) ! 2 or u) 2 


by integral ®. Similarly, the identical gain must come about for the magnetic 
field, integral ®, through Vr/> since no coefficients appear. 


G ( 6 , 0) > !H(x\ y', 7 1 ) I 2 or ! 70 2 



In the far-field, only O(IA) 1 terms remain, and the wave number k = 2n A 
establishes the fact that 


G(0, 0)> |H(x # , y', 7 ')! 2 t oj 2 


It is anticipated, therefore, that 



is invaria nt with frequency for cophased apertures of the type assumed, observed 
on-axis, providing a check on computer programs. Useful decisions, concerning 
the need for double precision, sampling etc., can be made by insisting on this o 
invariance of the integral. 
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When the cophased aperture assumption is removed, it can be seen that 


I ( n x II j ) <//dS 


leads to complex phasor jxdygons which are now not straight lines, but may 
have contributions tending to annihilate to some extent. They yield reduced field 
strengths, and imply gain degradation of the system. In the far-field region, each 
transverse component , E f of the scattered field leads to a phasor polygon. 

No attempt will be made here to relate the envelope or shape of the scattered 
radiation pattern to absolute signal strength consistent with energy conservation 
principles via the integral, over all space, as required by the definition of gain. 
The implications of the superposition principle (a linear concept) and the notion 
of energy conservation (a square-law for fields) are not peculiar to the Kirchhoff- 
Kottler formulation. 1 

Having discussed the frequency dependence of integrals ( 2 ) and ( 7 ) in the far- 
field for a cophased aperture system, integrals ® and ® are examined as 
these complete the formulation under the assumption of perfect conductivity. 

Since V0 implies the first power of o, it follows that the Kottler contour integral 
is frequency inde pendent due to the presence of the coefficient o ~ 1 outside of the 
integral. The surface charge integral (?), however is without any coefficient and 
varies directly as the first power of oj or frequency. To sum up: 


E( x ' , y * , z' ) 


I® (^°) + + !(4)( a ') 

invariant 


H(x'.y'.z') = Iq(o,) 


for simple cophased apertures observed on-axis. 


THE INTERMEDIATE NEAR-FIELD 

Returning to the general formulation, and assuming a perfectly conducting sur- 
face, only the integrals (T), ©, and @ of the electric field and the integral © 


‘Kef. 10, pp, 101- 10*»- 
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of the magnetic field remain and are examined here in some detail. The behavior 
of the sum of these integrals in the intermediate near-field region is of interest, 
as is the special case when the observer recedes to infinity. Also, prior to be- 
coming engrossed in computational aspects of the problem several questions 
might be asked concerning the Kirchhoff-Kottlcr formulation: 

(1) Does the solution for the scattered field satisfy Maxwell's equations? 

(2) Can one expect the solution to satisfy the radiation conditions at infinity? 
(i.e. Is there a net radial field component?) 

(3) Is the Kottlcr integral an appended "correction" term (approximation) 
based on physical intuition or argument? 

(4) What is the behaivor of the Kottler integral for closed surfaces since the 
contour C was originally for open surfaces ? 

The answers to these questions can be obtained without resorting to evaluation 
involving particular geometries, illumination, etc., even though the task is 
tedious and involves numerous manipulations on the field expressions. 

For the reader who is interested primarily in the application of the theory, and 
to a lesser degree in the mathematical arguments underlying it, the answers 
are tabulated here. The derivations can be found in the Appendix of this document. 

(1) The Kirchhoff-Kottler formulation leads to solutions £(x* , y' , z' ) f 
H( x' , y ' , z ' ) which satisfy Maxwell's equations. 


— r/ lj — U U 

V x E = - JjT ’ V x H = J + J7 ( free space) 


(2) The radiation conditions are satisfied by the Kirchhoff-Kottler formulation 
and the radial components of integrals ® , © , and (?) annihilate for the 
remote observer of the electric field. The magnetic field is inherently 
transverse at infinity since the radial field is precluded by a cross-product 
involving a radial vector. 


1 i m r E( x ' , y ' , z ' ) is finite, 

r -• ® 
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z') is finite, 


l r x E(x' , y' , z') -H(x' 
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The proof lies in a development 1 which shows that both E(x' , y ' , z ' ) 
and H(x' , y ' , z' ) are divergenceless, so that 


V' • E(x', y\ z') = V' • H(x\ y', z') = 0 


Further evidence is presented by a separate proof due to Sancer 2 who 
shows that 


p/ ¥ * v > . < \ — bn i.i y 

E( x , y , z ) - ^ 


If ^r) * I" (nxH l )e ,k> ds 

Js, 


as r 


00 


H(x' , 


y', O 



as r - 00 . 


The latter development substantiates in detail an assertion by Silver 3 
that the far-zone field is transverse by the current distribution method 
(i.e. by the Kirchhoff-Kottler method). 


Ref. 1, pp. 469-470; and Appendix A, this document 
“Ref. 7, p. 143; and Appendix B, this document 
3 Ref. 2, p. 149 
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(3) The Kottler integral, which is introduced via a physical argument that 
envisions a discontinuity in a sheet current is shown by Sancer to evolve 
naturally using the free-space dyadic Green's function. 


G - (i ♦^-VVg) 

where 


g = 


± 

4*7 


4 Tt 


e 


- j k 



It appears that the Kottler contour term is more than a "correction" 
term that incidentally leads to the satisfaction of the Maxwell equations. 


(4) The Kottler integral for a closed surface vanishes identically for arbi- 
trary illumination distribution because a dyad analogue to Stokes' 
theorem 1 is applicable and clearly shows that the contour itself vanishes. 

That is 




* ( H, V./,) 


) V0H, • di - 0 


1 


1 


See Appendix A for the vector and dyadic form of the curl theorem of Stokes on an open surface- 
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since 




1 



The intermediate near-field formulation is now examined. As shown in Fig. 1, 
an arbitrary reflecting surface, feed-point, and point of observation are con- 
sidered. The phase and magnitude of the wave front arriving at the surface is 
affected by a factor of the form e" * Vp, /p ' , and accounts for the path length 
traversed in passing from feed to reflector. This is already stated implicitly 
in the theory when E, and H, are specified at y x . 



Figure 1. The Intermediate Near-Field Region 


Only three paths are traced from the feed to the observer to suggest the process 
being described. It can be seen that the distance and direction from feed to 
surface must be assumed unique, and the same is true of the distance and direc- 
tion from surface to observer. A factor e“ jkr /r is evidently required for the 
latter part of the path. Over the entire path the product (e - 1 kp //c)(e" J kr /r) 
describes the amplitude and phase of a propagated wave . 

In the Kirchhoff-Kottler formulation the history of the fields emanating from the 
surface is of concern, and is given by the quantities 4> = e* J kr /r and V<// 
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- ( jk ♦ l/r)«/>l ( . These factors are under the integral sign, in general, since 
each differential area is uniquely affected. Since the charge distributions are 
not associated with propagation, but fields, V0 appears instead oi 0 above. By 
this device, the formulation introduces terms which descend as 0( 1/r )* and 
0(l/r) 2 and simultaneously achieves the (local) radial character of those fields 
arising from charge distribution. See integrals (T) and ( 4 ) . Only the quantity 0 
appears for integral (J), and the orientation of this field component is determined 
by the direction of a (local) unit tangent, (nx H^/Inx ll t s T. The (total) field 
given by integral © is purely transverse (in a local context) since, effectively, 
the tangent 7 is crossed into a radial vector 1 r contained in 

To conclude, it is noted that the backseattered fields of the Kirchhoff-Koltler 
formulations are formed in Cartesian components E x , E y , E, and later resolved 
into spherica* components E r , E^, E,. From Figure 1 is can be seen that sup- 
posed radial fields, which can only attract or repel along a line of action between 
a differential area ds and the observer, can have an E 0 and/or E^ resolute with 
respect to the moving spherical triad of basis vectors. Likewise, supposed trans- 
verse fields arising from a differential area in such a manner that they are 
orthogonal toVi/ (or l r ),can have radial resolutes relative to the spherical sys- 
tem. The 1 r , 1 < 9 , 1 . basis vectors are shown in Figure 1 to illustrate this fact. 
The notion that the 1„, 1 . vectors lead to a description of a transverse field 
is not particularly useful, therefore, in the intermediate near-field region. 
Mathematically the resolution into spherical components is correct but it does 
not appear to aid the understanding of the physical assets of the problem. Also, 
since fields other than those synchronous and orthogonal fields associated with a 
spherical wave are now obtained in the solution, it appears that dual reflector 
systems may be measurably affected by induction as well as propagation terms. 


SUMMARY 

This document assumed the Kirchhoff-Kottler diffraction formulation presented 
in the literature and discussed some of the implications of theory for different 
specializations. A consideration of the behavior of the formulation when the 
reflectors are idealized open or closed surfaces, and when the observer is at a 
large distance from the scattering object, for example, provided valuable connec- 
tives between the mathematical theory and the physics of the problem. Recovery 
of the intrinsic impedance of free space from the ratio of E( x' , y ' , z ' ) to 
H( x ' , y ' , z ' ) , and the gain dependence on v 2 , as well as the vanishing of the 
radial fields at infinity, not only satisfied intuitive ideas on diffraction, but also 
provided means for verifying a computer program. Satisfaction of Maxwell's 
equations by the general solutions for the scattered fields was also demonstrated. 
Finally, the derivation of the Kirchhoff-Kottler formulation in a direct manner 
via the dyadic Green's function afforded an interesting departure from the his- 
torical development. 
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In conclusion, the objective of this document was to provide an overview of the 
Kirchhoff-Kottler formulation and stimulate further Interest in diffraction studies 
by their method of analysis coupled with present-day computer technology. 

The modular Fortran diffraction program develo|)ed at Goddard Space Flight 
Center by the Antenna Systems Branch has been used in a large variety of ap- 
plications including monopulse configurations, thermally distorted satellite 
antennas, large arrays of source elements and dual-reflector systems. While the 
mathematics and theoretical discussions may at times appear abstract, the 
ultimate objective has been that of determining the extent of the Inherent soundness 
of the Kirchhoff-Kottler approach for simulating objects in the physical world. 

The combined effect of wave curvature, surface curvature and size, and observer 
range and angle on this 88-year old theory are still not generally known. 

"...to pursue mathematical analysis while at the same time turning one's back on 
its applications and on intuition is to condemn it to hopeless atrophy." (R. 

Courant) 1 Increased computer capability should erJiance applications and thereby 
clarify the range of validity of the theory. 
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APPENDIX A 


STRATTON'S PROOFS (FIELD DIVERGENCE AND 
MAXWELL'S EQUATIONS) 1 


Stratton has shown that the Kirchhoff-Kottler formulation for E( x ' , y ' , z ' ) , 
consisting of four integrals, has zero divergence and leads to waves which are 
transverse at great distances from the sources. He has also shown that 
V' xH(x\ y ' , z ' ) = (d/d t ) E( x ' , y ' , z ' ). It is stated that H( x ' , y ' , z ' ) has zero 
divergence, and V' x E( x ' , y ' , z ' ) = - (OB/d t ) ( x ' , y ' , z ' ). The former proofs 
are carried out on the general formulation, and the conclusions are therefore, 
also valid when the boundary conditions for a perfect reflector (cr -«) are 
invoked. 

In this appendix the proofs of Stratton are given in detail for both V' • E(x' , y ' , z ' ) 
and V' • H( x' , y ' . z ' ) . For completeness, V' x E( x', y', z') = - (d/dt ) B(x', y', z') 
and V' xH(x', y', z') ~(d/d t) E(x', y\ z') are both verified. Stratton's operator 
V v* is not used, and associations such as V • (H t are treated as the diver- 
gence of a dyad (Rj Vi/# ) , since neither a vector dot- or cross-product is in- 
dicated. All of the terms arising through the use of vector identities are written 
out in detail. The arguments which subsequently vanish some of these terms are 
then presented. In some instances identities could not be found for operators, 
vectors, and dyads and only identities for vectors and dyads were available. In 
such instances detailed expansions were actually written as indicated, but will not 
be given here as they are easily verified. It was not assumed that an operator 
could be substituted for a vector in these cases. Some of the more tedious 
expansions involved as many as 72 terms to prove an identity relationship, and 
these are not reproduced here. Even so, the derivations still tend to be tedious; 
however, they provide valuable insight cerning the behavior of the Kirchhoff- 
Kottler formulation. 


TABULATION OF IDENTITIES 


& 

III 

V</> • A + <£V • 

(I) 

/*N 

>1 

X 

col 

B'Vx^-A'VxB 

(II) 


‘Ref. 1, pp. 469-470. 
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V x (V</>) * 0 


(III) 


I 


(V X A) 


n (IS 



d? 


curl theorem of Stokes 


V x ( </>A ) = V x A + V</> x A 


( A x B) • C = A • (BxC) 


Vx(AxB) = (AV • - A • V) B - (BV • - B • V) A 


di • D 



V x Dds 


dyad analogue of Stokes’ 

theorem 


A • (BxD) = (AxB) • D where D is a dyad 


L 


V 


A dv 


= £ b . 


A ds 


divergence theorem of Gauss 


• D dv 



n • D ds 


dyad analogue of Gauss' theorem 


(IV) 

(V) 

(VI) 

(VII) 

(Vin) 

(IX) 

(X) 

pci) 
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THE TRANSVERSE ELECTRIC FIELD 

Repeating the general formulation for the scattered electric field: 


E(x , y', z') = r— — 


Ilf 

joj€ 4n j 1 


dZ 


4n J xH,) 0 f (n x Ej x V// + ( n • Ej V0 J ds 


The transverse character of the field at great distances follows from 1 


V' 


E(x\ y\ z') 





H, * dZ) • V0 f (H, 


<ll)V‘ • (V0) 


(I) 


r - 

j 0 • (nx H,) + 0V' • (nxH,)lds (I) 

*' *i 


" 4^ J [^0 • V' x (nx Ej - (nx E,) • V' x (V0)] ds 

Sl (II), (III) 

“ J [ V ' ( n ‘ Ej * (V0) + (n • E,) V' • (V0)J ds (I) 

^ The Roman numeral designates the tabulated identity used. 
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Since V' - -V when the operators apply to p or its derivatives ; 1,2 and since the 
illumination distribution is not affected by a change of coordinates on the part 
of the observer ( x ' , y ' , z ' ), 


V' • E(x\ y\ z' ) 


= - 4 - 4- 1 

JO)€ 477 J 


V 2 0 H, * dl 


[ [W(n xH ,) 1 Vp ♦ (n * E,) V 2 </>] ds 
J S . 


±1 

ja>€ 



0 H, * dt 


kl 

4 77 



n • E,) pds 


+ 



• Vi// ds 


] Ref. 1, p. 469, p. 169. 

^Stratton uses <p = e + * ltr /r and Silver uses p - e ’’“/r for the solution to the homogeneous wave 
equation p + k 2 p - 0. In this document </> is reserved for a spherical coordinate or a generic 
scaler with identities. Take p - e above. 
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• cK 


I (0VxHj + V0 x H , ) * nds (IV), (V) 

J"! 


- (V0xH t ) • n 


(n x Hj) * V// 


(VI) 


V' -E(x\ y ' , z') = 


= - "£7 izj W, -d-t --^rj(n • E,) *ds 


r , - 7 r 

- -4^- 4* '/'H, • d-t t-^rj 


+ ~ ± 4 ^ ~ | 0V X H , • n ds 


Then 


V' • E( x ' , y', z') = 0 


since 


V x Hj 





- jcje Ej e~ iwt . 
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THE TRANSVERSE MAGNETIC FIELD 

Repeating the general formulation for the scattered magnetic field: 


H( x , y\ z') = 


joj/n 4 .t j 




In J [j"«( n kEJ 0 - (n x Hj ) x V0 - ( n • Hj ) V0 jds 


The transverse character of the field at great distances follows from 


V' • H(x\ y\ z') = 


1^71^ v '(Ei • d^) • 70 + (Ej • dt)7' 


W' (I) 


477 j* jo:e jv* 0 • ( n X Ej) 4 07* ■ j n x E, ) Jds 


(I) 


1_ 

4tt 


[70- V' x(n x Hj j 

•'S. L 


-(nxH,)-7'x(V,/,)]ds<II), (III) 


1_ 

4^7 


| f" V' ( n • H, ) • V0 + ( n 
Jsj 




Hj) 74.* (V0)Jds 


(I) 
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H'x\ y\ z') = ' ?M E t •«*) 


J x E i) * W' “(n • Hj )V 2 «/>Jds 


k 2 


a( 

•' c 


k 2 

■^E -M --^ 


j (n*H,)^di 

Js. 


¥ [ (n»E,) 

Js. 


V<// ds 


[v x ^ 0 Ej)J ’n ds = j) (//Ej’ d't = J (</ / VxE 1 + V 0 x Ej * n ds (IV), (V) 


“ x E^ • n = InxEj) • V 0 


(VI) 


H( x , y\ z') = 


k 2 

jOft 


4tt| ^ E 1 


TP k 2 


d-£ “ 


4?f 


[ ("•».) 

Js. 


0 ds 


jo;6 r 

-4F J ^ 


d^ - 


OJ€ 


4?r 


j 0 (VxEj) • n 

•J S - 


ds 
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Then 


V* • H(x\ y\ z') = 0 


since 


V x E, 


j9 

at 


( B i e ' i Jt ) 


jcfiHj 


e-jo* . 


THE CURL OF THE SCATTERED MAGNETIC FIELD 
Repeating the general formulation for the scattered magnetic field: 


H(x\ y\ z') = 


1 

jaft 


l 


WE, 


M 


I 

( 




^-(nx H,) x Vi/i-|n • Hj) 



ds 


Maxwell’s equation VxH = J + (a D/at) (free space) is satisfied as follows: 


V' xH(x\ y\ z') = [(E, -d^) V'x V^+V' (E, -d^)xV 0 ] 


(V) 


^ | x(nxE,) +(7'0)x(nx EjjdstV) 
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An 



n ,) * ♦ ( n x h t ) 



(VII) 





Wp + V' ^ n • Hj) x V/-jds(V) 


where 


V' x^JnxH,) x VV/J 


(n<Hj) V' • V0 


- ^ n x HjJ • V' J V0 - V0 V' *(exHj) + • V') (nxH,) 


k 3 (n xH,) i/< + | nx Hj) • Wi// . 


(VII) 


It can be verified that 


| ( n x Hj ) * W \p ds 

*'*1 



n • V x Hj) ^ds 



n x V • ( Hj V0)ds 


by treating ( H, V <p) as a dyad 1 and writing out the 72 terms of the equality. 
Derivatives of the illumination distribution with respect to the operator y do not 
vanish. 


! Ref. 6, Appendix 3; Ref. 11, Chapters V, VI, and VII; Ref. 12, p. 420, p. 520; Ref. 13, p. 46, p. 113; 
Ref. 14, See Index. 
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Using the dyad analogue to Stokes' theorem, 




ft • (H, W/) 


| n • V x (Hj Vt//j ds 



n < V * ( Hj V// )ds 


(VIII) 


The last equality' was verified, term-by-term, by expanding ( FI t Vv) , and also 
follows from identity IX when B - V. In general, operators cannot replace vectors 
or dyads in identity relationships. Therefore, 



• Wi p ds 


“ )Cl>€ 



n • F.j) V// ds 



V0H, 


• d-6 


and 


V' x H(x\ y', z‘) 


■ 


Vi/zHj • dl 


| [( n x E i) x W' + j xH,) /+ (tr E,) w] 


ds 
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so that 


V' x H( x ' , y ' , z ' ) = - jo>f E( x ' , y ' , z ' ) 


and Maxwell’s first equation is satisfied by the Kirchhoff-Kottler formulation. 

THE CURL OF THE SCATTERED ELECTRIC FIELD 
Repeating the general formulation for the electric field: 


E( x ' , y\ z' ) 


J_ JL f 

l 4 J 


W'H, 


cl7 


l_ 

+ rr 



Hj »/> + (nx Ej) x ^ + 




ds 


Maxwell’s equation V x E = - ) is satisfied as follows: 


E( x 1 


O 



(Hj • dl) V' x V0+ V' (H, -d^jxV^J 


(V) 


4^ 




(V' *p) x ( n x Hj) J ds 


(V) 
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where 



V' x ^ n x Ej) x Vi//J = ( n x Ej ) V* • V</> n x E 4 ) • V'J 


V<// 


“ V// V' • Inx Ejj (V// 


= k 2 |n xEj /< + (iix E,) • VV0 


It can be verified that 



W0 ds 



xEj ] V// ds 



n x V 


* Ej) x V//J ds(V) 


V') (n x E,) 


(VII) 


• (E Vi//) ds 


by treating (Ej Vi//) as a dyad and writing out the individual terms of the equality. 



Using the dyad analogue to Stokes' theorem, 


| W/fEj • d 4) = (d4-E,)W/ = | d* • (E, V0) 


| n • V x [Ej Vi//j ds = I n x V • 
■'•l •'•i 


(E t w)ds.(vni) 


Therefore, 



Wt// ds 





n • H, 



V0 E, 


• d-f 


and 



|ja)e^n xEj ~ (n x Hj) x Vi/> - • Hj j V<//j ds 

1 



so that 


t 

V' x E( y', z 1 ) = joj/i H(x', y', z') 


and Maxwell's second equation is satisfied by the Kirchhoff-Kottler formulation. 
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APPENDIX B 


SANCER'S DERIVATIONS (THE FREE-SPACE 
DYADIC GREEN’S FUNCTION ) 1 


Sancer has shown that the complete Kirchhoff-Kottler formulation forL(x' y' z' ) 
can be obtained using the free-space dyadic Green’s function, thereby avoiding 
separate physical or intuitive arguments to obtain the contour integral. Sancer 
has also manipulated the general formulation in such a way as to establish the 
fact that the radial fields vanish at infinity as required . 2 The special case of a 
closed surface is treated in terms of two open surfaces using the vector form 
of Stokes’ theorem. 3 

In this appendix, the derivations of Sancer are not reproduced since that author 
uses conventional notation, operators, etc. throughout. This appendix augments 
Sancer’s article through a discussion of certain dyadic identities, arbitrary 
vectors, and operators in spherical coordinates. References to texts containing 
special identities are given, and dyad analogues to expressions known from 
vector analysis are identified. Wherever the expressions are too involved for 
complete presentation, the method of obtaining the desired result is simply 
outlined. 

The dyadic Green’s function 


G 




K 


where 


g 



exp. ( ik r 



^ef. ?. pp. 141-144. 

^Ref. 2, p. 149 and p. 88. Silver states that "The effect of the boundary-line distribution is 
therefore to cancel the longitudinal field component introduced by the surface charge and current 
distribu:ions. M Reference is then made to a volume integral argument for the far- zone fields. The 
Kottler integral is not explicity involved. 

3 Ref. 15, p. 352. 
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is fundamental to Sancer’s development. 1 One of the first relationships involving 
G is 2 


(LG) • V' = L(G-V') 


where V'(7' ) is an arbitrary vector field, and the operator L - Vx Vx - k 2 . The 
given relationship is verified in a straightforward manner by expanding both 
sides, beginning with the curl of the curl of the dyad G . 

Now V x G 



~ * 
i fly + 




i j 


+ G 


X Z 


rs a 

ik 


y x 


j i + G yy j j 


+ G 


y * 


jk 


kl + G *y kj 


A A 1 


/ f)G * 


1 ^ G x 

3G A A, 

( dG y 

<*G\ 

\dy 

<9z / 

' + j ("a? ■ 

' «9x / + k 

( ax ‘ 

' a y / 


where 

k 

•A 

k 

k . 


G = G i + G i + G 

x XX xy J xz 


G = G i + G j + G 

y yx yy J y z 


G = G i + G j + G 

z zx zy J zz 


1 

2 


Ref. 6, p. 221 Van Bladel makes extensive use of the dyadic Green’s function in solving a larg 
variety o f mechanics and electromagnetics problems. 

Ref. 7, p. 142, Equation (2.8). 


e 
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Then V x (VxC) 





/aG v 

£S,\ 

_d_ 

/rVG x dG\l 

•l?i 

f \~dt 

dy ) 

dz 

V dz ‘ dx/J 



fL) 

3 (S 


[d z \ d y 

■ c)z! 

0X \ dx 

" a y /J 


"_a_ f^x 



r dG 

Z 

a6 ,)l 

_(9 X \ d Z 

lx / 

' dy' 

K dy 

dz/_ 


The analogy with vector analysis is obvious. In the expression 1 for V x Vx V , 
replace vector V by the dyad G . 

The left-hand sid§ of the equation can now be expanded, and the result dotted 
intoV' to form (LG) • V'. The right-hand side of the equation is then treated 
by forming the vector 


G • 


V' 



+ G 


xy 



+ G 


X Z 



+ 



V 


t 

X 


*Ref. 7, p. 490, identity 59. 
instead of d"\^/dxdz> 


+ G V'+G V')T+(g V'+G V'+G V '1 

yy y y z z / J \ z x x zy y zz z J 

A correction is required. The last factor should read d" /dydz 
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7T> 


and writing out VxVx(G-V'). The identity 


(LG) • V' - L(G • V) 


is verified by this process, however, it appears that the arbitrary vector field 
V' (r ) is a constant vector field relative to the unprimed partials. Writing only 
the i component of the equality here, to save space, the left-hand and right-hand 
sides of the relationship yield 


2 G 0 2 G 2 G 

yjc xx xx 


V x \dy Ox ’ d z 2 


Oy 2 


o 2 c \ 

*x ) 

{ Oz Ox ) 


+ V ' 


(- 


«? 2 G 0 2 G 
yy *2 


y\?ySx dz 2 


d 2 G H 2 G 

U , DL 

Oz Ox 


) 


//) 2 r <9 2 G <9 2 G d 2 G \ 

V ' _ G y i - 11 li + LL 

1 \ dydx dz 2 dy 2 OzOxi 


Oy Ox 


(G V ' 
\ y * * 


G V ' 
y y y 


G V') 
y* * / 



v *' 


+ G. 


xy 


V 


+ G 



dz 


r G V ' 

. XX X 


G V ' 

xy y 


+ G 


X z 



+ 


0 2 

01 H 



+ G 


* y 
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It ts noted that the dot product of a dyad with a vector Is non-commutatlve in 
general. The curl of the curl of a dyad is a dyad, and a dyad dotted with a 
vector is a vector. 

After applying the divergence theorem of Gauss,' Sancer writes 


(. 


V' 



which expresses the fact that G is symmetric. As expected, an expansion of the 
equality shows that the symmetry refers to the principal diagonal symmetry of 
the matrix representation of G, so thatG ( = G , , when the i, j are the row- 
column designation. Sancer next employs a relationship 


V x (G • V' ) = Vg x V' 


which is easily verfied by expanding the left-hand side and writing out the expres- 
sion G - (I + l/k 2 VV)g. Since 


_ rr o* 2 ~~ d 2 d 2 

4 ‘i dx dy ' *k<9xdz 


~ ~ d 2 d 


rl 2 A 2 

A A 1/ A A (/ 


dy dx jj 2 + ^ dy d z + ^ * dz dx + ^ j d z dy 4 


kk 


dz 2 ' 


Vx(G-V) = i [^r (gv, ) -JT (*v y )] 


i[s7 (s v „) ' ^T(8 V .)] + k '[£ («v y )-£(gv„)] • 


1 


See the identities X and XI in Appendix A. 
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I 


Then the identity is verified when V' is an arbitrary constant vector field rel- 
ative to the unprimed partlals. 

After a considerable amount of manipulation , 1 which is straightforward, but 
tedious to verify, Sancer arrives at a result 2 that is written in his notation as 


F, . ) 




+ (n' x H' ) + — V'V' g • (n x H)\ ds 

CD 


The integrals CD and (T) of Stratton’s formulation can be identified. Further 
manipulation which employs various artifices such as adding and subtracting 
the same internal, finally results in 3 


E(r) 


E i (r) + j , {( n x E k' ) * V 'fc 

^ s i 

CD 


+ ioftg (n' x H k ') + • E k ')Jds,' 

(D CD 




^The triple scalar product and identity IX of Appendix B are used. 
^Ref. ?, p. 142, Equation (2.17). 

^Ref. 7, p. 143, Equation (2.24). 
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which can be identified as the original formulation of Stratton.' The existence 
of the contour integral is es|)ecially significant since it was obtained directly via 
the free-space dyadic Green’s function. The intermediate expression, above, 
will be used subsequently to investigate the radial fields in the far -zone. 

Returning to an argument by Sancer 1 concerning the contour integral for closed 
surfaces. 


v' -f-Z7 


V R ( Hj ” H 2 ) <U 


is said to vanish when H is continuous on the surface due to Stokes' theorem for 
a closed surface. The argument 's subtle since, if H is continuous, H = M, = H 2 = 0 
;*md it is easy to lose sight of the fact that it is the contour which vanished. 3 

The following alternative viewpoint is presented to avoid the a priori statement 
that the surface integral vanishes, since it is the vanishing of the contour integral 
which establishes the far-less obvious vanishing of the surface integral. Ob- 
viously, the presence of the arbitiary field vector V allows Sancer to employ 
the vect or form of Stokes' theorem. The behavior of the contour integral can 
be determined without the use of V' by using the dya d analogue to Stokes' 
theorem. 

From Stokes, in Stratton's notation, 


n • V x(Hj Vi// ) ds , 


visualizing a simple open surface Sj initially. The contour is by common con- 
vention, oriented so that the right-hand rule holds for the normal to the surface. 
If a second simple open surface should exist such that the boundary c, is also 


£ 


V<// H j • d'C 


■{ 


d-t • (H, ^ ) = 


I 


! Rcf. 1, p. 469. 

“Ref. 7, p. 143, Equation (2.23). 
3 Ref. 15, p. 352. 
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the boundary - c 2 , it follows that 


1 


n • V x ( Hj Vi// J ds 


1 * * 2 * * 



d' • (H, 


> d*, • (H, W) 

'l 


because the surfaces and contours are orientablo. That isd-ij = - cU 2 . Then 
the concluf ion is that the integral over the closed surface also vanishes, 1 
necessarily, since th^ contour has been annihilated. 

Sancer’s intermediate result 


A <\i 2 - (h, V0) = A (di l ♦ d>e 2 ) (H, 7/,) = 0 

J c 2 " c 


E(r ) 


M r > 




E')x7*g 


CD 


+ ioqj.% (n' x H' ) 

CD 



(n x 


H')| 


ds 


leads to the conclusion that the Kirchhoff-Kottler formulation exhibits transverse 
fields only in the far-zone limit. Recalling the dyadic Green’s function, 



^Ref K, p. 338. Several examples of orientable and non-orientaLie surfaces are given with sketches, 
and dir«*ct reference to Stokes’ theorem is also made when S = Sj + S 2 as above. 
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Sancor writes 


E(T) 


E, (r) ♦ 


ICi +1C 


» k r 


4rrr 




(n X H>‘ 



ike ' kr _ f _ _ / 

ds' ♦ — - — n x I (n* x E' )e~ iy * r * r ds' , 
4rrr r J » 


after Injecting the far-field approximation e” ikn r ' r , and rewriting G in soherical 
coordinates in the limit as r -oc . 

The GF( r > follows from the spherical operator 


7 


f _£ - 1 ^ - 1 ± 

* r Sr + Iff r 57? + 1<^ r sin# Si ’ 


and the basis vectors of the moving spherical triad 


A A - 

l r = sin ' cos 0 i + sin# sin# j + cos # k 


\q - cos ^ cos ^ i + cos# sin# j -sin# 1 



s in # i + cos # j 
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Then 


W R 



+ 


d ( ~ d r \ ^ d ( ~ _1 <9 r \ >9 / ~ 1 <9 r \ 

r sin 5 d<$> r <9 r/ + r s i n 0 <9^ rdf 1 ) 4 r s in (? <9$ \*<*> r s in 0 d<t>) 


* ~ <9 2 R 

= 1 1 — 

r ' dr 2 


\e & g *<*> ** / g • 

r % r r dr 



i r i - >< 2 ) 


+ 1* 1 


e l e 


/ i k 1\ ~ ~ / i k 1 \ 

^ r r 2 )^ ^ ^ \ r ~r 2 ) R 


From this representation of Wg, 1 



1 A — 

The various authors use different symbols to represent the unit radial vector. Here i f = a f (Sancer) 
(Van Bladel) ~r (Newton) which is usually clear from the context. Also, g = e ,kr /r (Sancer) and 
g = e _jk 7r (Van Bladel). 
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In the limit 1 as r - 00 , 


0 


Cf (r) = 


fi ‘ n , “,)« 


as used by Sancer. When G ( r ) operates in a dot product against an integral of 
the form 


I- 


H' )e 


• i ku r 


dS' 


the result can be regarded as the entire vector field minus the radial part of 
that field via (I - a f a r )g. Alternatively, the result can be regarded as the dot 
product of a dyad G® ( r ) which has only \ e and 11. components in a matrix, 
with a field of 1 r , \ c ! \ , values. This is ecpjivalent to taking only the transverse 

components of the field. 

To sum up, Sancer has shown that the Kirchhoff- Kottler formulation has only 
transverse fields in the far-zone limit since (1) 

(n* x E' ) x V g ds 



is necessarily transverse (but vanishes anyway when the conductor is perfectly 
reflecting), and (2) the radial part of the integral for electric sheet current is 
annihilated by an integral which is equivalent to the boundary curve integral on 
electric charge plus the integral on surface-charge distribution. 


‘Ref 7, p. 221 and Ref. 17, p, »*"» 
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